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Abstract. Motivated by recent experiments where interference patterns behind 
a grating are obtained by accumulating single photon events, here we provide 
an electromagnetic energy flow-line description to explain the emergence of such 
patterns. We find and discuss an analogy between the equation describing these 
energy flow lines and the equation of Bohmian trajectories used to describe the 
motion of massive particles. 
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1. Introduction 

The possibility of performing quantum interference experiments with low-intensity 
beams (i.e., one per one particle) of photons [1-3] and material particles [4-6] has 
intensified the theoretical search of the topology of the photon paths [7-9] and particle 
trajectories [10-14] that describe the process behind the interference grating. The aim 
of all the proposed approaches is to simulate the appearance of the interference pattern 
by accumulation of single-particle events. 

In Bohmian mechanics one may simulate this process for material particles. 
Bohmian trajectories follow the streamlines associated with the quantum-mechanical 
probability current density and, therefore, reproduce exactly the quantum-mechanical 
particle space distribution in both the near and the far fields [10-12]. Alternatively, 
the emergence of the interference pattern in the far field has also been simulated by 
sets of rectilinear trajectories characterized by the momentum distribution associated 
with the particle wave function [13,14]. In the far field, the distribution of momentum 
components along Bohmian trajectories is consistent with this distribution [14]. 

In this paper, we show how to determine electromagnetic energy (EME) flow lines 
behind an interference grating, where the components of the magnetic and electric 
vector fields satisfy Maxwell's equations. These fields are expressed in terms of a 
function which explicitly takes into account the boundary conditions imposed by the 



Electromagnetic energy flow lines as possible paths of photons 



2 



grating. The EME flow lines are then determined after numerically solving the path 
equation arising from the Poynting energy flow vector. In particular, here we show 
EME flow lines behind gratings consisting of different number of slits. These sets 
of lines supplement those presented by Prosser [7] for both a semi-infinite plane 
and gratings with one and two slits. The EME flow lines show that the energy 
redistribution behind the grating until reaching the Fraunhofer regime. In particular, 
it is of interest the process that corresponds to multiple-slit gratings, where we 
can observe the smooth transition from a Talbot pattern in the near field to the 
characteristic Fraunhofer peaks in the far field. 

It is tempting to conclude from the results obtained that the motion of an eventual 
photon wave packet thus represents an energy flow along a group of flow lines. This 
conclusion is supported also by the fact that the path equation for the EME flow 
lines has the same form as the equation for the quantum flow associated with material 
particles. This explains why there is a complete similarity in interference phenomena 
with photons and material particles. Experimentally, the final interference patterns 
as well as the processes of their emergence are analogous [1-6]. 



2. The complex Poynting vector and the equation for the EME flow lines 

The diffraction of electromagnetic radiation by a grating is described by the solution 
of Maxwell's equations in vacuum that satisfy the grating boundary conditions [15]. 
We consider the simplest solutions of Maxwell's equations: harmonic electromagnetic 
waves 

H(r,i)=H(r)e--*, (1) 
E(r,i) =E(r)e- <w *. (2) 

The physical electric and magnetic fields are obtained by taking the real parts of the 
corresponding complex quantities. The space-dependent parts of these fields (which 
are complex amplitudes) satisfy the time-independent Maxwell equations 

V x H(r) = -iwe E(r), (3) 

V x E(r) = iw/i H(r), (4) 

V-H(r) =0, (5) 

V-E(r) =0. (6) 

From these equations it follows that both fields, E(r) and H(r), satisfy the Hclmholtz 
equation 

V 2 H(r)+fc 2 H(r) = 0, (7) 
V 2 E(r) + fc 2 E(r) = 0, (8) 

where k = u/c. 

The EME flow lines are now determined from the energy flux vector, which is the 
time-averaged flux of energy, given by the real part of the complex Poynting vector [16] 

S(r) = lRe[E(r)xH*(r)]. (9) 

Note that, since the flow of energy goes in the direction of the Poynting vector, the 
EME flow lines can then be determined from the parametric differential equation 
dr 1 S(r) 



ds cU (r) ' 



(10) 
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were s is a certain arc-length along the corresponding path and U(r) is the time- 
averaged electromagnetic energy density 

U(r) = \ [e E(r) • E*(r) + W H(r) • H*(r)]. (11) 

3. The Poynting vector in the case of two dimensional diffraction by a 
plane grating 

To describe simple diffraction experiments, we will consider that the grating is on 
the XZ plane at y = 0, with the incident plane harmonic wave traveling along the 
^-direction. Moreover, we assume the problem to be completely independent of the 
z-coordinate (i.e., very long slits along this coordinate). In order to encompass all 
possible cases of polarization, we express the magnetic and electric fields before the 
grating as a superposition of two waves: H polarized, for which the magnetic field is 
along z (A components), and E polarized, for which the electric field is along z (B 
components) [17]. That is, 

H(r, t) = U(r)e- iut = [Ae lky e z + Be^e^eJ e~ iut , (12) 

E(r, t) = E(r)e- 1 "* = — \-Ae iky e x + Be^e lky e z ] e~ tuJt . (13) 
e uj 

Here, <$> is the phase shift between the x-and z-component of the field, and the constants 
A and B are real. For (j> = or tt, the incident wave is linearly polarized, while the 
cases A = B, with <fr = ±7r/2, describe circular polarization. The cases A ^ B, with 
4> ^ 0,7r and ±7r/2, describe elliptic polarization. 

As shown by Born and Wolf [17], the i?-polarization and _ff-polarization 
components satisfy two independent sets of equations since the problem is independent 
of the z-coordinate. This implies that the solution that corresponds to an incident 
wave, given by (fT2"j) and (fT3"| . diffracted by a grating can be expressed as 

H = ik^Be^ ^e x - ik^Be^ ^-e v + Aibe„ (14) 
dy dx v Y ' y ' 

iA dif> iA dif> kB id> , 

E= ^ e, -^e y + e*fye„ 15 

eoiv ay cqui ox cqlo 

where tp(x,y) is a solution of the Hclmholtz equation 

VV + fcV = 0, (16) 

which satisfies the grating boundary conditions. From (|14[) and (|15[) . we can now 
express the components of the time-averaged Poynting vector in terms of the constants 
A, B and </> and the function ip(x,y), as 

s * = ^+*A* 9 i -*•£). <") 

12 , o2^ (.,.&!>* _ L *dll> 



*=^ A ' + ^[*w- r w> (18) 

S - 1 ABzin^t 9 ^ 9 ^* d ^* d ^\ ( 19 ) 
z 2eowfc \ dx dy dx dx ) 

In the case of a linearly polarized initial wave ((f) — or tt), it follows from 
equation (fT9|) that the z component of the Poynting vector vanishes. Thus, the EME 
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flow lines will remain confined on the XY plane and, from the parametric differential 
equation (|10[) . we find that the differential equation 



dy _ S y _\ dy dy 



dx S x ( dtp* „ dtp 
dx dx 



(20) 



will determine the eventual photon paths, which are obtained by numerical integration. 
As can be noticed, the topology of the EME flow lines in the case of linear polarization 
is independent of the constants A and B, and therefore independent of the direction 
of polarization. 

On the other hand, for circular and elliptic polarizations the z-component of the 
Poynting vector is non zero, this leading to important differences in the properties of 
the corresponding EME flow lines, which will not be planar, as infers from the non 
vanishing z-componcnt in equation (|10p . The properties of EME flow lines in these 
cases are beyond the scope of the present work and will be explored in a forthcoming 
one. 



4. Flow lines behind a specific grating 

In order to plot flow lines for a specific grating we need explicit expressions for the 
magnetic and electric field behind a grating. Traditionally, explicit solutions have 
been written using the solution of the Helmholtz equation in the form of the Fresnel- 
Kirchhoff integral [17]. By making the appropriate approximations, the solution was 
then transformed into the expressions valid in the Fresnel and Fraunhofer regions, 
respectively. Talbot effect and Talbot-Laue effect were also explained [18] by making 
the appropriate approximations and transformations of the Fresnel-Kirchhoff integral. 

If the x component of the wave vector satisfies the relation k 3> k x , the solution 
of the Helmholtz equation may also be expressed as a superposition of transverse 
modes (STM) of the field multiplied by an exponential function of the longitudinal 
coordinate [19]. As shown by Arsenovic et al [20], this form is equivalent to the 
Fresnel-Kirchhoff integral. For an incident plane wave falling on the grating at y = 
and lying on the XZ plane, the STM form of the solution reads as 

p iky roc 

1>(x,y) = Z— c(k x )e* k **e-* k *y/ 2k dk x , (21) 



27T J- 

where the function c(k x ) is determined by the incident wave and the grating boundary 
conditions as 

1 f°° 

c(k x ) = ^= iP{x,0)e- ik * x dk x . (22) 



Within the approximation k 3> k x , one finds 



dtp 
dx 



dtp 
dy 



— « iktp, 
dy 



(23) 



and, therefore, the EME density (fTTj) becomes proportional to \tp(x,y)\ 2 , i.e. 



U(r) = -n (A 2 + B 2 )\tP(x,y)\ 2 , 
which is obtained from ([Ti"]). (fl~5|) and (fTT]) . 



(24) 
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Figure 1. EME flow lines behind a Ronchi grating (<5 = d/2) with N = 5 slits, 
A = 500 nm, d = 20A = 10 /im and L T = d 2 /\ = 200 jim. 



Assuming that the grating is completely transparent inside the slits and 
completely absorbing outside them, we have that ip(x,0) — for points outside the 
slit aperture and ip(x,0) = vp in (x,0) for the inside points, where tp in (x,0) is the field 
component of the incident wave. In the case of an incident plane wave propagating 
along the y-axis, tpi n (x,0) is a constant. For a grating with N apertures separated a 
distance d and all with the same width <5, a simple integration [19] renders 

J_ IT sm(k x 6/2) am(Nk x d/2) 
Cl xj "V2^ViV k x S/2 sin(M/2) - 1 j 



In figures [T] and [21 we have represented the EME flow lines behind Ronchi gratings 
with 5 and 30 slits, respectively. The unit along the y-axis is the so-called Talbot 
distance, Lt = d 2 /A, which gives the repetition period behind a grating of the 
diffracted wave [12]; the unit along the axis parallel to the grating (x-axis) is the 
period d this grating. As seen, the topology displayed by the EME flow lines is very 
similar to that displayed by Bohmian trajectories for massive particles [10-12]. Note 
that the quantum density current 

J = — (W#* - #*W) (26) 
2m 

determines the particle trajectories in Bohmian mechanics through the equation of 
motion 

Tt - i ^ 

dt p 

where p — W*4 r . The latter equation is analogous to equation (|10p . though this 
analogy is not totally complete: in the case of particles with a mass m, the relation 
y = hkt/m holds, whereas the same is not true for photons. When the wave function 
^ does not depend on z, J has only x and y components, 

(») 



2m \ dx dx 

\J/ V]/* 

2m \ dy dy 



J y = — i^-tr- ■ ( 29 ) 




Figure 2. EME flow lines behind a Ronchi grating (6 = d/2) with N = 30 slits, 
A = 500 nm, d = 20A = 10 fim and L T = d 2 /X = 200 fMm. (a) EME lines from 
the upper part of the grating (y > 0). (b) EME lines from six central slits of the 
grating. 



Equations (|28[) and (f2"5)) have the same form as equations (|17p and (HHJ). Thus, 
the EME flow lines in the case of a linearly polarized incident wave are similar to 
the Bohmian trajectories of massive particles described by a two-dimensional wave 
function. 

5. Emergence of the interference pattern by accumulation of single 
photon arrivals 

In the case of massive particle the modulus square of the wave function describes 
the distribution of particles at a distance y from the grating in the far field after 
(theoretically) an infinite number of particles has reached the detector (at y). This 
theoretical result has been nicely confirmed by a new generation of experiments which 
use low-intensity beams of particles. In these experiments, the final interference 
patterns are built-up after particles accumulate gradually one by one at a scanning 
screen [2,3,5,6]. Numerical simulation of particle arrivals, assuming that they 
move along de Broglie-Bohm's trajectories [10-12] and MD trajectories [13] describe 
theoretically this process. This means that a trajectory-based interpretation completes 
the standard interpretation of the wave function, where a picture in terms of single 
events is missing. 
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(a) 








Figure 3. Histogram of the number of trajectories ending at various points along 
the ce-axis at a distance y = A.3Lx for four different values of the total number of 
photons: (a) 100, (b) 1000, (c) 2000 and (d) 5000. Here diffraction is produced 
by a two-slit grating and the initial conditions (positions along the two slits) for 
the photons data are chosen at random. The red line is the plot of the function 
\ip(x, y)\ 2 . It is seen that at the chosen distance y = 4.3Lt, maximum of the 
distribution at x/d = 8.6, associated with k x d = An, is absent. This follows 
from (25) because for N = 2, d = 25, the second interference maximum, being at 
k x d = An coincides with the first diffraction minimum. 



Analogously, one can proceed in the same way with photons assuming that they 
move along the EME flow lines described in the previous section. This is illustrated in 
figure[3j where we have plotted histograms (blue dots) obtained from the accumulation 
of photons for a two-slit diffraction experiment (with A = 500 nm, d = 20A and 
S = d/2). In particular, the histograms have been made considering equally spaced 
bins (with a width of 0.9<i) at an observation distance y b — A.ZLt from the plane of the 
grating, where the Fraunhofer pattern is already well converged (this happens when 
the observation distance is further away than the so-called Rayleigh distance [21], 
which in this case is yn ~ (25 + d/2) 2 /4nX ~ 0.2Lt <C y b)- As can be seen, as 
we move from panel (a) to (d) the histogram data approach better and better the 
smooth (red) line which represents the EME density U(x, y b) given by equation (f24f 
as the number of photons per time unit increases, as also seen in the experiment [2,3]. 
It is interesting to note that the photons (or, equivalently, the initial positions of 
their paths) distribute randomly along the distance covered by each slit aperture and, 
therefore, their arrival positions at y = A.ZLt will also be random. However, they will 
accumulate in accordance to U(r) because of the guidance condition given by equation 
(I10|) . which can be alternatively expressed [17] as 

S(r) = U(r)v, (30) 

where v is a sort of effective vectorial velocity field which transports the EME density 
through space in the form of the EME density current. The vector field v is always 
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oriented in the direction of the wave vector k. Thus, before the grating it is aligned 
along the y direction and behind the grating will depend on the particular point (x, y) 
where the field is evaluated, becoming almost constant along some specific direction 
only within the Fraunhofer regime. 
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